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求解 并 行 机 拖 期 与 能 耗 成 本 优化 调度 的 混合 教 一 学 算法 “ 
王 永 琦 ， 吴 发， 江 潇潇 ， 王 春 姐 
(上 海 工程 技术 大 学 电子 电气 工程 学 院 ， 上 海 201620) 


摘 要 : 针对 加 工时 间 可 控 的 并 行 机 调度 ， 提 出 了 一 类 考虑 拖 期 与 能 耗 成 本 优化 的 调度 问题 。 首 先 对 调度 问题 进行 了 
问题 描述 , 并 建立 了 整数 线性 规划 模型 以 便于 CPLEX 求解 。 为 了 快速 获得 问题 的 满意 解 , 提出 了 一 种 混合 教 -学 算法 。 
结合 问题 的 性 质 ， 设计 了 编码 与 解码 方法 以 克服 标准 教 -学 算法 无 法 直接 适用 于 离散 问题 的 缺点 。 同 时 ， 构建 了 基于 变 
邻 域 搜索 的 局 部 搜索 算 子 以 强化 混合 算法 的 搜索 性 能 。 最 后 ， 对 加 工时 间 可 控 的 并 行 机 调度 问题 进行 了 仿真 实验 ， 测 
试 结 果 验 证 了 本 文 构建 的 整数 线性 规划 模型 和 混合 算法 的 可 行 性 和 有 效 性 。 
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Hybrid teaching-learning-based optimization algorithm for optimizing tardiness and energy 
cost on parallel machine scheduling 


Wang Yongqi, Wu Fei, Jiang Xiaoxiao, Wang Chunyuan 
(School of Electronic & Electrical Engineering, Shanghai University of Engineering, Shanghai 201620, China) 


Abstract: For the parallel machine scheduling with controllable processing times (PMS-CPT) , this paper proposed a Scheduling 
problem for optimizing tardiness and energy cost. First, this article developed an integer linear programming (ILP) model for 
the proposed scheduling problem in order to facilitate the CPLEX solver. Second, this paper established a hybrid teaching- 
learning-based optimization (HTLBO) algorithm to quickly obtain satisfactory solutions. This algorithm adopted a new coding 
and decoding method according to the nature of the proposed parallel machine scheduling problem, which overcame the 
weakness of standard teaching-learning-based optimization algorithm that could not be directly applied to discrete problems. 
Meanwhile, a local optimizer was proposed based on variable neighborhood search (VNS) , which aimed to enhance the 
performance of the hybrid algorithm. Finally, the simulations were conducted to solve instances of PMS-CPT. The experimental 
results demonstrate the feasibility and validity of the proposed ILP model and hybrid algorithm . 
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optimization algorithm 


人 的 重要 方面 PH。 因 此 ,进行 并 行 机 拖 期 与 能 耗 成 本 优化 调度 方 
折 的 研究 具有 重要 的 理论 意义 和 工程 应 用 价值 。 
作为 当前 实际 生产 制造 过 程 中 的 一 类 典型 的 调度 问题 ， 并 PMS 问题 已 被 证 明 为 NP-hard 问题 ， 快 速 、 高 效 的 算法 设 
行 机 调度 (parallel machine scheduling, PMS) 问 题 具有 相当 广泛 。” 计 已 经 成 为 这 一 领域 的 重要 课题 。 相 关 的 求解 算法 可 以 归 


的 应 用 背景 , 包括 云 计 算 、 半 导体 加 工 、 汽 车 制造 等 行业 [ 
统 的 PMS 问题 致力 于 优化 与 作业 完成 时 间 相关 的 指标 ， 包 括 。” ”包括 动态 规划 、 分 枝 定 界 等 算法 , 这 类 算法 能 够 获得 PMS 


最 大 完工 时 间 、 提 前 期 / 拖 期 等 。 随 着 客户 定制 化 的 快速 发 展 以 ”的 准确 解 ， 但 随 着 问题 规模 的 扩大 ， 这 类 算法 所 需 的 计算 时 间 


。 传 以 下 三 类 : 精确 算法 、 启 发 式 算法 和 智能 算法 四。 精确 求解 算法 


纳 为 


T 


问题 


及 全 球 能 源 价格 的 攀升 ， 准 时 化 与 节能 调度 成 为 当前 车 间 调 度 。 将 急剧 增长 。 启 发 式 算法 利用 调度 规则 来 实现 加 工 任务 在 不 同 


t 
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录用 稿 


机 器 上 的 分 配 及 排序 ， 能 够 快速 获得 PMS 问题 的 近似 最 优 解 ， 


但 这 类 求解 方法 对 问题 设置 具有 较 严格 的 限制 。 近 年 来 ， 随 着 
智能 算法 的 快速 发 展 ， 学 者 们 已 经 将 各 类 智能 算法 运用 到 相关 
各 
提 


PMS 问题 的 求解 中 来 ， 包括 遗传 算法 、 禁 忌 搜 索 算法 、 粒 子 群 
算法 和 模拟 退火 算法 等 ， 这 类 算法 能 够 在 满足 时 间 要 求 的 前 


下 获得 PMS 问题 的 满意 解 四 。 然 而 ， 上 述 智能 算法 存在 控制 参 

数 多 难以 选取 、 收 和 敛 速度 慢 ,迭代 后 期 易 陷 入 局 部 最 优等 缺点 。 

为 此 , 本 文 构建 了 混合 教 -学 算法 以 快速 有 效 地 求解 PMS 问题 。 
教 - 学 优化 算法 (teaching-learningbased optimization, TLBO) 


是 Rao 等 学 者 提出 的 一 种 新 型 的 群体 智能 优化 算法 站。 该 算法 


通过 模拟 实际 师 生 授课 过 程 中 的 “教学 ”和 “学 习 ” 两 个 过 程 
来 实现 问题 的 优化 ， 以 上 两 个 过 程 在 TLBO 算法 中 分 别 定义 为 


“教学 阶段 ”和 “学 习 阶 段 ”"”。TLBO 算法 具有 控制 参数 少 、 


收 


敛 速度 快 等 优点 ， 并 在 复杂 函数 优化 、 神 经 网 络 训练 、 资 源 项 


算法 应 用 于 并 行 机 调度 问题 的 研究 。 


目 调 度 等 多 个 方面 都 获得 了 成 功 应 用 -1 同时 , 鲜 有 文献 将 该 


基于 以 上 研究 ， 本 文 针对 加 工时 间 可 控 的 并 行 机 调度 ， 提 


出 了 一 类 考虑 拖 期 与 能 耗 成 本 优化 的 调度 问题 。 构 建 了 该 问题 
的 整数 线性 规划 模型 以 利于 CPLEX 求解 ， 并 提出 了 混合 教 一 


学 求解 算法 以 快速 获得 中 大 规模 问题 的 满意 解 。 最 后 进行 了 仿 


法 的 可 行 性 和 有 效 性 。 


1 ”问题 模型 


1.1 问题 描述 


真 测试 ， 力 求 验 证 本 文 构建 的 整数 线性 规划 模型 和 混合 优化 算 
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王 永 琦 ， 等 : 求解 并 行 机 拖 期 与 能 耗 成 本 优化 调度 的 混合 教 一 学 算 ; 


tj :机 器 ;以 速度 5, 加工 订单 j 所 需 的 时 间 , 即 1 = 0,/ 5 
决策 变量 : 
xu :表示 0-1 变量 , 当 订 单 j 在 机 器 i 上 位 置 1 处 以 速度 8， 


加 工时 , 取 值 为 1; 否 则 ,为 0; 


bp, :机 器 i 上 位 置 1 处 的 订单 开始 加 工 的 时 间 。 


mn nh 大 


Min FD YD 0) 
/=! 


i=1 j=1 1=1 s=1 


m 


a O) 


i=l 1=1 s=1! 


nn 


>》 Dx S11zsigsm, 1<1<n (3) 
j=1 3=1 


b= Dt xi lSigm, 1<1<n-l (0 


j=1 s=1 


b, 


il+tl 


了 >b,, + fis* Xi —d—xis)M —d; 
l<i<m,l<s<u, 1<j<n,l<l<n (5) 


T,>0, 1<l<n (0) 


Xi E{0,1}, 1<igm,; 
l<s<u, 1<l<nls jn (7) 
其 中 ， 式 (1D) 定义 了 目标 函数 ， 该 目标 函数 包含 两 项 成 本 : 所 有 
订单 的 拖 期 惩罚 成 本 总 和 ; 所 有 订单 的 加 工 能 耗 成 本 总 和 。 约 
束 式 (2) 表 述 每 一 个 订单 只 能 安排 在 一 台 机 器 的 一 处 位 置 上 以 一 
个 确定 的 速度 进行 加 工 ， 约束 式 (3) 表 示 每 台 机 器 上 的 一 个 位 置 
处 至 多 安排 一 个 订单 的 加 工作 业 ; 约束 式 4) 计 算 了 机 器 ;2 上 位 


假设 有 m 台 机 器 负责 n 个 订单 的 加 工作 业 , 每 个 订 利 


了 各 


7 处 订单 开始 加 工 的 时 间 ; 约束 式 (5)(6) 用 于 计算 订单 j 的 拖 


自 的 交 货 期 和 单位 时 间 拖 期 惩罚 成 本 ， 每 台 机 器 有 多 种 不 同 的 
加 工 速度 ， 不 同 加 工 速度 对 应 的 单位 时 间 的 能 耗 成 本 不 同 。 
本 文 致力 于 优化 并 行 机 调度 问题 的 拖 期 与 能 耗 成 本 总 和 ， 


需要 决策 以 下 三 方面 的 因素 :不 同 机 器 所 负责 加 工 的 订单 集合 ; 
各 个 机 器 上 订单 的 加 工 顺序 ; 各 个 机 器 上 所 有 订单 的 加 工 速 度 。 


1.2 ”数学 模型 


PMS-CPT 问题 的 整数 线性 规划 模型 。 
m :机 器 总 数 ,机 器 编号 以 ; 表示 ; 
n :订单 总 数 ,订单 编号 以 j 表示 ; 


其 于 以 上 问题 描述 ， 在 此 定义 以 下 符号 以 构建 本 文 提 


出 的 


1 :表示 机 器 上 的 加 工 位 置 编号 ,各 个 机 器 均 设置 n 个 加 工 


位 置 , 即 1<7<m 
u :表示 机 器 ;的 加 工 速度 总 数 ,速度 编号 以 s 表示 ; 
4 :订单 j 的 工作 负荷; 
dj :订单 j 的 交 货 期 ; 
Ci :订单 j 的 完工 时 间 ; 
7 :订单 j 的 拖 期 , 即 T = max{0,C) 一 a,} 
oj :订单 j 单位 时 间 的 拖 期 成 本 ; 
5S; :机 器 ;第 s 个 加 工 速度 的 取 值 ; 


。 :机 器 ; 以 速度 * 运转 时 产生 的 单位 时 间 的 能 耗 成 本 ; 


期 7 ; 式 (7) 给 出 了 自 变 量 的 取 值 范围 。 
2 ”标准 教 一 学 优化 算法 

TLBO 算法 利用 种 群 来 表示 问题 的 解 ， 当 前 种 群 中 最 优秀 
的 个 体 称 为 教师 ， 其 余 个 体 称 为 学 生 。 算 法 的 邻 域 搜索 部 分 包 
括 “ 教 学 阶段 ” 和“ 学习 阶 段 ” 两 个 部 分 。 在 教学 阶段 ，TLBO 
算法 利用 教师 个 体 和 群体 均值 个 体 对 学 生 个 体 进行 教学 ， 进 而 
实现 种 群 的 进化 ， 在 学 习 阶 段 ，TLBO 算法 从 当前 群体 中 随机 
地 选择 两 个 个 体 ， 并 采用 当前 个 体 向 较 优 个 体 进行 学 习 的 策略 
实现 种 群 的 更 新 。 标 准 TLBO 算法 的 实现 流程 归纳 如 下 : 

a) 初 始 化 TLBO 算法 的 控制 参数 和 待 优化 问题 的 参数 。 包 
括 种 群 规模 学 生 数 目 p ， 问 题 变 量 的 维 数 p ,各 个 维度 自 变量 
的 取 值 范围 [x/,xs] 。 
b) 构 建 初始 种 群 ， 并 计算 每 个 个 体 的 适应 度 值 ， 个 体 的 初 


台 化 公式 为 x =X +rand(0,)) (到 - 芭 ) 。 其 中 ， yand(0,1) 表示 


[0, 了 上 的 随机 数 。 


©) 教 学 阶段 。 首 先 ， 计 算 当 前 群体 所 有 个 体 每 一 维 的 均值 
元， 进而 得 到 平均 个 体 ， 元 = 伺 , 元 ， 台 }。 其次， 对 于 当前 种 


人 


群 中 的 任 一 个 体 x* ， 利 用 个 体 更 新 公式 得 到 新 个 体 w 。 最 后 ， 
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录用 稿 
比较 新 旧 个 体 xw 和 x 的 适应 度 值 ， 并 采取 贪 禁 保 留 策略 ， 即 保 


留 适 应 度 值 较 优 的 个 体 ， 放 弃 适 应 度 值 较 劣 的 个 体 。 个 体 公 式 


为 : X=x+rand(0,D) (x -0-x)， 其 中 * 为 教师 个 体 ， 即 当前 


种 群 中 的 最 优 解 , 教学 因子 8 表述 随机 整数 , 取 值 集合 为 {1,2} 。 

d) 学 习 阶 段 。 首 先 ， 随 机 从 经 历 过 教学 阶段 的 种 群 中 选择 
两 个 个 体 , 分 别 记 做 xt 和 x? ;其 次 , 选择 x! 和 x? 中 较 优 秀 的 个 
自 定 待 优化 的 问题 为 最 


体 进 行 学 习 。 以 f(.) 表示 适应 度 函 数 ， 
小 化 问题 ， 个 体 更 新 公式 为 
Se 若 f (x') < f (x) 
Xx+rand(0,D) (x 一 x) 否则 
e) 算 法 终止 条 件 判 断 。 若 满足 算法 终止 条 件 ( 达 到 最 大 迭代 
次 数 或 规定 的 计算 时 间 等 ), 则 终止 TLBO 算法 , 并 输出 当前 最 
优 解 ， 否 则 ， 转 ¢)。 


3 ”改进 教 -学 算法 


编码 与 解码 
标准 TLBO 算法 采用 实数 编码 方式 ， 无 法 直接 适用 于 离散 
问题 ， 故 需要 结合 PMS-CPT 问题 的 性 质 在 编码 与 解码 方面 进 
行 重新 的 定义 和 设计 。 

已 知 问题 参数 : 机 器 总 数 m ， 订 单 总 数 n 和 每 个 机 器 的 加 
工 速 度 总 数 w (1<i<m), 则 本 文 提出 的 并 行 机 调度 问题 可 采用 
1xn 维 数 组 进行 编码 表示 ， 各 维度 编码 数值 的 取 值 范围 为 
[m+ 了]) ,编码 的 第 j 维 对 应 订单 ) 。 第 j 维 编码 的 整数 部 分 对 
应 订单 了 所 分 配 的 机 器 编号 ; 第 j 维 编码 小 数 部 分 第 


3.1 


Cu 


1~|logw*wm2…m | 用 于 确定 订单 j 的 加 工 速度 ， 第 


10 


[log 2 ml 位 至 最 后 用 于 确定 订单 j 在 机 器 上 的 加 工 


顺序 。 
解码 过 程 : 首先 ， 对 编码 进行 整数 取 整 ， 以 此 确定 各 个 订 
单 对 应 的 机 器 编号 。 其 次 ， 取 各 维 编码 小 数 部 分 的 第 


构建 各 个 机 器 的 轮 盘 赌 概率 模型 以 


1~ | 10g lial | 位 ， 


选择 各 订单 的 加 工 速度 ; 针对 机 器 ;, 已 知 该 机 器 的 速度 类 型 总 
数 为 w ， 则 当 相 应 编码 数值 处 于 区 间 [(s-D/w,s/w) 时， 订单 
的 加 工 速度 取 值 为 5, (1<i<m，1<s<u)。 最后, 依据 分 配 到 


3 


同一 机 器 上 的 订单 对 应 的 编码 数值 的 第 [iog'wr2…m]Hl 位 


至 最 后 的 数值 大 小 ， 进 行 升序 排列 以 确定 订单 加 工 顺序 。 
图 1 给 出 了 编码 与 解码 示意 图 ， 算 例 参 数 为 : 机 器 总 数 


m=2 ， 订单 总 数 n= 每 个 机 器 的 加 工 速度 种 类 u,=3 ( l1<i<m )。 


因此 ， 编 码 的 各 个 维度 的 取 值 范围 为 [3) ， 编 码 的 整数 部 分 表 
示 订 单 -机 器 分 配 情况 , 编码 的 小 数 部 分 第 一 位 用 于 确定 订单 加 
工 速度 ， 小 数 部 分 第 二 位 至 最 后 用 于 确定 订单 加 工 顺序 。 已 知 
初始 编码 {2.629, 2.812, 1.254, 2.827, 2.265, 1.195} 。 解 码 结果 归 


纳 如 下 : 
订单 1 在 机 器 2 上 以 速度 5,, 进行 加 工 ; 
订单 2 在 机 器 2 上 以 速度 9,; 进行 加 工 ; 
订单 3 在 机 器 1 上 以 速度 5 进行 加 工 ; 
订单 4 在 机 器 2 上 以 速度 9,, 进行 加 工 ; 
订单 5 在 机 器 1 上 以 速度 9 进行 加 工 ; 
机 器 1 的 订单 加 工 顺序 为 {3,5}; 
机 器 2 的 订单 加 工 顺序 为 {2,4,1}。 
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图 1 编码 与 解码 示意 图 


3.2 基于 VNS 的 局 部 搜索 算 子 


与 其 他 群 智能 优化 算法 相 类 似 ,标准 TLBO 算法 存在 着 迭 


代 后 
建 了 基于 3 
部 搜索 算 子 ， 力 求 


期 易 陷入 局 部 最 优 、 寻 优 性 能 不 足 等 缺点 。 为 由 
> 邻 域 搜索 (Variable neighborhood search，VNS) 的 局 
曾 强 标准 TLBO 算法 的 搜索 性 能 。 


上 ， 本 文 构 


VNS 算法 通过 系统 地 改变 邻 域 结 构 集 以 拓展 当前 算法 的 


搜索 范围 ， 进 而 获得 局 部 最 优 解 ; 


同时 ， 基 于 此 局 部 最 优 解 ， 


再 次 系统 地 改变 邻 域 结构 集 来 进 


步 拓展 搜索 的 范 


居 ， 力 求 找 


到 另 一 个 局 部 最 优 解 (4。 下 图 给 出 了 VNS 算法 搜索 过 程 的 示 


意图 。 
当前 最 小 值 
Ax) 对 应 邻 域 
x \ 


局 部 最 优 解 
对 应 邻 域 


图 2 VNS 算法 示意 图 


给 定 解 x={x,%,…,Xxp}，VNS 算法 的 核心 在 了 


@ 局 部 最 优 解 。” 率 全 局 最 优 解 


确定 x 的 邻 


域 结 构 和 搜索 流程 。 给 定 循环 次 数 及， 本 文 构建 的 
局 部 搜索 算 子 可 表述 为 : 


< 


于 VNS 的 
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a) 令 天 =1， 转 b)。 

b) 令 z=1， 转 c)。 

oO 车 qz<D， 转 由 ; 否则 ， 转 全 。 

gd) 令 X=+sign.rand(0,]D).B.m， 转 6)。 

9 令 d dtl， 转 0)。 

人 计算 wr 的 适应 度 值 fx) ; 若 f(x) < FoD， 则 令 r< mw， 
转 g)。 

多 令 1 和 +1， 若 人 < 及， 转 b); 否则 ， 结 束 VNS 搜索 。 
其 中 ，d 定 义 了 邻 域 变换 公式 ， sign 表述 随机 整数 ， 取 值 集合 
为 {-1,1} ，B 为 邻 域 控制 参数 取 值 为 (0,1] ，m 为 模型 参数 ( 即 机 


3.3 ”改进 教 -学 算法 流程 
结合 以 上 描述 ， 给 出 求解 PMS-CPT 问题 的 HTLBO 算法 


流程 ， 归 纳 如 下 : 
a) 初 始 化 算法 的 各 个 参数 ， 包 括 种 群 规 模 、 算 法 终止 条 件 


稚 


b) 采 用 随机 方法 构建 初始 种 群 ， 并 计算 适应 度 值 。 

9 计算 当前 种 群 的 平均 个 体 , 并 利用 “教学 阶段 ”的 邻 域 更 
新 公式 更 新 种 群 。 

gd 利 用 “学 习 阶段 ”的 邻 域 更 新 公式 更 新 种 群 。 

e 针 对 “教学 阶段 ” 和 “学 习 阶 段 ” 个 体 更 新 失败 的 解 ， 采 
用 基于 VNS 的 局 部 搜索 算 子 进行 个 体 的 二 次 更 新 操作 。 
人 D 若 算法 满足 终止 条 件 ， 则 终止 迭代 搜索 并 输出 当前 最 优 
解 ， 否 则 ， 转 co)。 


4 ”仿真 实验 及 分 析 


4.1 算 例 设计 
于 难以 从 现 有 的 文献 中 获得 与 本 文 提出 的 并 行 机 调度 问 


i 


naXiv:201805.00379v1 


i 


ch 


题 类 似 的 测试 数据 集 , 因此 , 本 文 依据 Hall 和 Posner 学 者 提出 
的 机 器 调度 问题 测试 数据 生成 方法 构建 了 本 文 的 测试 算 例 D3。 
本 文 构建 的 测试 算 例 分 为 : 小 规模 算 例 和 中 大 规模 算 例 ， 
# 计 21 个。 其中, 小 规模 测试 算 例 能 够 通过 CPLEX 进行 求解 ， 
将 所 得 的 最 优 解 用 于 混合 教 -学 算法 的 性 能 分 析 ; 中 大 规模 算 
例 则 无 法 在 合理 的 计算 时 间 内 借助 CPLEX 获得 精确 解 ， 因 此 
采用 其 他 同类 型 智能 算法 所 得 的 满意 解 进行 混合 教 -学 算法 的 
性 能 分 析 。 对 于 小 规模 算 例 ， 机 器 数目 me {3,4,5} ， 订 单数 目 
ne{5,10,15} ; 对 于 中 大 规模 算 例 ， 机 器 数目 me {3,4,5}， 订单 
数目 ne {30,50,80,100}。 对 于 算 例 (m,n)， 相关 参 数 取 值 如 下 : 
引 订 单 j 的 工作 负荷 4 在 均匀 分 布 U(,100) 上 产生 ; 订单 
j 单位 时 间 的 拖 期 成 本 系数 w 在 均匀 分 布 U(0,1) 上 产生 ; 
b) 各 机 器 加 工 速度 类 型 总 数 w=5， 机 器 ;第 s 个 加 工 速度 
的 取 值 8, 在 均匀 分 布 (5,10) 上 产生 
0o) 机 器 i 以 速度 s 运转 时 产生 的 单位 时 间 的 能 耗 成 本 e 在 
均匀 分 布 (0,1) 上 产生 ， 且 加 工 速度 越 低 ， 对 应 的 e, 值 越 大 ; 
订单 j 的 交 货 期 d 在 均匀 分 布 Uy(0.5.D,2.D) 上 产生 ， 
其 中 5 的 表达 式 为 
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在 MATLAB2016a 平台 上 进行 编程 测试 ， 计 算 机 参数 为 主 
频 1.6 GHz、 内 存 8 GB、Intel Corei5-8250UCPU，PMS-CPT 问 
题 的 整数 线性 规划 通过 CPLEX-MATLAB 接口 借助 CPLEX 
12.4 进行 求解 。 针 对 各 个 测试 算 例 ，HTLBO 算法 进行 30 次 独 
立 仿真 测试 。 在 小 规模 算 例 中 ，HTLBO 算法 的 种 群 规模 设 为 
30， 最 大 计算 时 间 设 为 10s， 参 数 天 取 值 为 5， 参数 1 取 值 为 
0.2; 在 中 大 规模 算 例 中 ，HTLBO 算法 的 种 群 规模 设 为 80， 最 
大 计算 时 间 设 为 1308， 参 数 到 取 值 为 135， 参数 8 取 值 为 0.1。 
4.2 测试 结果 
首先 ， 利 用 CPLEX 和 HTLBO 算法 对 小 规模 测试 算 例 进 
行 了 求解 ， 测 试 结果 见 表 1。 其 中 ， p* 表示 CPLEX 所 得 的 精 
确 解 ， 后 一 列 为 CPLEX 求解 各 个 算 例 所 需 的 计算 时 间 ; 石和 
五, 分别 表示 HTLBO 算法 在 30 次 独立 仿真 测试 中 所 获得 的 最 
优 解 和 最 劣 解 ， 忆 表示 获得 最 优 解 的 实验 次 数 占 总 实验 次 数 
〈 即 30 次 ) 的 百分比 。 
由 表 1 的 测试 结果 可 知 , CPLEX 能 够 利用 本 文 构建 的 整数 
线性 规划 模型 获得 PMS-CPT 调度 问题 的 精确 解 ( p* )， 据 此 验 
证 了 整数 线性 规划 模型 的 有 效 性 。 但 CPLEX 所 消耗 的 计算 时 
间 较 大 ， 算 例 (5,15) 达 到 了 2815s。 对 于 各 个 小 规模 测试 算 例 
(m,n) ;本 文 构建 的 HTLBO 能 够 获得 问题 的 准确 解 ( 即 =F* )， 
且 30 次 独立 实验 所 得 准确 解 的 百分比 ( PB) 不 低 于 83.74%; 同 
时 ， 最 劣 解 ( 己 ) 与 准确 解数 值 较为 接近 。 考 虑 到 HTLBO 的 计 
算 时 间 仅 为 10s， 远 小 于 CPLEX 的 计算 时 间 ， 据 此 可 验证 
HTLBO 算法 在 解决 小 规模 算 例 时 的 有 效 性 。 


表 1 小 规模 算 例 测试 结果 


CPLEX HTLBO 

(m,n) 

FF 时 间 /s 于 Ey BI% 
(3,5) 20.46 233 20.46 21.05 94.68 
(3,10) 38.52 412 38.52 39.44 93.16 
(3,15) 111.45 1686 111.45 113.21 91.16 
(4,5) 10.67 351 10.67 10.93 95.01 
(4,10) 28.14 1514 28.14 28.94 91.53 
(4,15) 87.52 1793 87.52 88.72 88.15 
(5,5) 13.35 564 13.35 13.80 91.80 
(5,10) 38.94 2382 38.94 39.44 88.09 
(5,15) 63.03 2815 63.03 63.69 83.74 


进一步 ， 进 行 了 中 大 规模 算 例 的 测试 。 考 虑 到 CPLEX 无 
法 在 合理 的 时 间 内 获得 这 类 测试 算 例 的 准确 解 , 因此 将 HTLBO 
算法 与 TLBO 算法 和 粒子 群 算法 (particle swarm optimization， 

PSO) 算 法 0 进行 比较 。 针 对 每 个 算 例 ， 各 算法 运行 30 次 ， 表 
2 统计 了 寻 优 结果 的 平均 值 ( 去) 与 标准 差 (c )。 此外, 变量 & 表 


录用 稿 


示 相 对 平均 偏差 百分比 , 针对 各 个 测试 算 例 (m,n), 令 三 个 算法 
共计 90 测 仿真 所 得 的 最 优 解 为 pr ， 变 量 & 的 计算 公式 为 
(F-F")/F'x100%。 
表 2 的 测试 结果 可 知 ， 融 合 了 基于 VNS 的 局 部 搜索 算 
子 的 HTLBO 算法 在 均值 ( F ) 和 相对 平均 偏差 百分比 ( x) 两 个 
指标 方面 均 要 优 于 TLBO 算法 和 PSO 算法 ， 据 此 可 知 HTLBO 
算法 的 搜索 结果 总 体 较 其 他 两 种 算法 更 优秀 。 同时, HTLBO 算 
法 30 次 独立 仿真 测试 结果 的 标准 差 (o 值 ) 更 小 ,这 说 明 HTLBO 
算法 具有 较 强 的 鲁 棒 性 ， 表 现 更 为 稳定 。 


5 ”结束 语 


瑟 


本 文 针对 加 工时 间 可 控 的 并 行 机 调度 ， 研 究 了 一 类 考虑 拖 
期 与 能 耗 成 本 优化 的 调度 问题 。 构 建 了 该 问题 的 整数 线性 规划 
模型 ， 并 借助 CPLEX 对 小 规模 问题 进行 精确 求解 。 同 时, 设计 
了 混合 教 - 学 算法 以 快速 获得 中 大 规模 问题 的 满意 解 。 对 21 个 
测试 算法 进行 了 仿真 测试 ， 实 验 结果 表明 CPLEX 和 混合 教 一 
学 算法 能 够 获得 小 规模 问题 的 精确 解 ， 且 混合 教 一 学 算法 的 计 
算 耗 时 较 短 ; 同时 相 较 于 标准 TLBO 算法 和 PSO 算法 , 混合 教 
一 学 算法 在 求解 中 大 规模 测试 算 例 时 表现 出 更 佳 的 性 能 ， 且 算 
法 鲁 棒 性 更 强 。 


| 


表 2 中 大 规模 算 例 测试 结果 
HTLBO 


TLBO PSO 


(m,n) 
F oO Kk/% F 他 Kk/% F wy Kk/% 


(3,30) 337.58 2.67 3.49 340.93 3.58 4.51 338.77 3.75 3.85 
(3,50) 812.57 7.66 4.21 819.78 12.50 5.13 817.76 10.32 4.87 
(3,80) 2354.43 24.10 5.35 2369.27 37.40 6.01 2360.39 32.28 5.61 
(3,100) 3138.57 34.17 6.30 3157.63 56.06 6.95 3149.36 43.80 6.67 
(4,30) 318.13 2.91 3.62 321.23 5.48 4.63 320.17 4.98 4.29 
(4,50) 625.26 7.61 4.02 632.52 12.18 5.23 629.51 10.23 4.73 
(4,80) 1679.45 12.09 5.09 1686.75 24.49 5.55 1683.92 20.81 5.37 
(4,100) 2183.21 16.46 6.66 2197.24 30.08 7.35 2192.75 34.46 7.13 
(5,30) 223.67 1.89 3.94 226.73 2.97 5.37 225.17 2.69 4.64 
(5,50) 450.66 3.89 4.47 457.76 7.16 6.12 454.09 6.93 5.27 
(5,80) 1142.95 12.50 5.79 1152.99 24.22 6.72 1145.40 20.08 6.01 


(5,100) 2208.56 20.45 6.72 2225.40 34.38 7.75 2217.73 31.27 7.62 
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